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ABSTRACT 


Gas  turbine  engine  design  is  based  on  simplifying  assumptions,  including  axi- 
symmetric  flow.  These  assumptions  break  down  in  the  real  machines  especially  when 
operating  near  to  stall.  The  reduction  in  operating  surge  margin  associated  with  current 
trends  to  develop  high  performance  military  fighter  aircraft  further  highlights  the  need  to 
develop  a  thorough  understanding  of  the  flow  in  rotor  blade  passages  and  especially 
when  operating  close  to  stall.  This  research  investigates  the  behavior  of  a  Transonic 
Compressor  Rig  (TCR)  as  it  approaches  stall  in  the  rotor-only  and  full-stage 
configuration  using  combinations  of  frequency  and  time  domain  analysis  in  the  subsonic 
(70%  rated  speed),  sonic  (80%  rated  speed)  and  transonic  operating  ranges  (90%,  95%, 
and  100%  rated  speed).  A  steady-state  analysis  of  the  pressure  measurements  across  a 
rotor  passage  of  the  TCR  was  conducted  to  identify  potential  pre-stall  indicators,  the 
structure  of  the  flow  field  within  the  blade  passages,  and  to  correlate  changes  in  flow 
structure  with  changes  in  pressure  measurements.  Further  investigation  using  flow 
visualization  techniques  within  the  blade  passages  demonstrated  how  the  TCR  rotor 
departed  from  an  idealized  rotor  when  the  axi-symmetric  assumption  holds  when 
operating  at  steady-state  and  near  stall  conditions.  This  identification  of  potential  stall 
indicators  and  characterization  of  the  time  variant  flow  field  within  the  passages  will 
impact  future  compressor  design  to  enable  operation  closer  to  stall  and  increased 
performance  across  an  operating  spectrum. 
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I. 


INTRODUCTION 


Rotational  stall  and  axi-symmetric  stall  and  compressor  surge  in  modem  aircraft 
gas  turbines  can  result  in  engine  efficiency  losses,  total  loss  of  thrust,  and  mechanical 
failure  due  to  high  cyclical  loading  of  the  relatively  thin  compressor  blades.  Both  stall 
modes  are  commonly  referred  to  as  “pop  stall”  because  of  the  flame  ejected  from  the 
engine  exhaust,  Figure  1.  This  is  often  caused  by  blade  wear  and  damage,  intake  of  high- 
temperature  gases,  and  extreme  flight  maneuvers.  Carrier-based  aircraft  are  particularly 
susceptible  to  this  type  of  engine  failure  because  of  their  extreme  operational 
requirements  and  the  likelihood  of  steam  ingestion  during  launch  from  the  carrier’s  steam 
catapult  system.  Because  of  the  reduced  stall  margin  associated  with  the  wear  and  age  of 
the  Navy’s  current  fleet  of  fighter  aircraft,  and  the  lack  of  engine  redundancy  introduced 
by  any  future  single  engine  aircraft,  a  thorough  understanding  of  compressor  stall  and  its 
effects  is  critical  to  mitigate  the  risks  associated  with  takeoff  and  operation. 


Figure  1.  Steam  Induced  Stall  of  F-18D  during  steam  catapult  assisted  takeoff  at 

NAWC  Patuxent  River 


Current  research  completed  at  the  Turbopropulsion  Laboratory  (TPL)  at  the  Naval 

Postgraduate  School  (NPS)  has  focused  on  the  steady-state  analysis  of  steam-induced 
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stall  of  a  transonic  compressor  rotor.  The  stage  (rotor  and  stator)  components 
implemented  in  this  study  were  designed  specifically  for  use  in  the  TPL  transonic 
compressor  rig  by  Sanger  at  the  NASA  Glenn  Research  Center  [1],  Rogers  [2],  Villescas 
[3],  Brunner  [4],  Payne  [5],  Zarro  [6],  Gannon  et  al.  [7]  induced  compressor  stall  by 
injecting  low  pressure  steam  through  the  rotor,  without  stator,  through  the  full  power 
range  (70%  -  100%)  The  compressor  performance  was  then  compared  to  pure  throttle 
induced  stall  to  identify  a  reduction  in  stall  margin  throughout  the  given  power  band. 
Koessler  [8]  and  Gannon  et  al.  [9]  furthered  this  research  using  high-pressure  steam  and 
identified  stall  precursors,  their  modal  oscillations,  and  an  additional  stall  margin 
reduction.  These  works  resulted  in  a  complete  performance  map  in  the  subsonic  (70% 
rated  speed)  and  transonic  ranges  (90%,  95%,  and  100%  rated  speed).  Gannon  et  al. 
further  extended  this  research  into  an  analysis  of  the  transient  stall  effects  of  the  rotor  in 
the  full  power  range  [10].  Using  these  previously  identified  performance  characteristics  a 
mapping  and  analysis  of  the  performance  of  the  full  stage  (rotor  and  stator)  was 
accomplished  by  Hurley  [11],  A  reduction  in  stall  margin,  when  compared  to  the  rotor- 
only  configuration,  was  identified.  This  indicated  that  the  stator  is  possibly  more 
susceptible  to  stall  than  the  rotor. 

This  paper  continues  the  full  stage  analysis  into  the  steady  state  and  transient 
phenomenon  present  in  the  subsonic  (70%  rated  speed),  transonic  (80%  rated  speed)  and 
supersonic  operating  ranges  (90%,  95%,  and  100%  rated  speed).  Additional  analysis  of 
the  test  rig  in  its  rotor  only  configuration  was  also  accomplished  using  an  updated 
instrumentation  array  and  procedure  when  compared  to  previous  rotor-only  research. 
Using  combinations  of  frequency  and  time-based  analysis  an  attempt  to  correlate  results 
across  methods  in  an  effort  to  validate  assumptions  and  interpretation  of  results  is  also 
presented. 
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II.  COMPRESSORS  STALL  THEORY 


Critical  to  a  compressor’s  ability  to  affect  an  increase  in  pressure  is  the  angle  with 
which  the  fluid  makes  relative  to  the  camberline  of  the  rotor  and  stator  blades.  Any 
angular  departure  that  the  fluid  makes  with  respect  to  the  tangent  of  the  camberline  is 
defined  as  the  incidence  angle,  Figure  2. 

Suction  Side 


Figure  2.  Angle  of  Incidence 

As  the  incidence  angle  increases  areas  of  low  pressure  can  form  on  the  suction 
side  of  the  blade,  Figure  2.  The  rapid  acceleration  of  the  flow  around  the  leading  edge 
must,  at  some  point,  begin  to  decelerate.  During  this  deceleration,  the  boundary  layer 
will  begin  to  grow  and  may  eventually  separate  resulting  in  a  region  of  reversed  flow. 

Since  a  boundary  layer  often  acts  like  a  physical  boundary,  the  sudden  increase  in 
layer  thickness  creates  an  effective  decrease  in  the  passage  area  between  the  blades.  This 
area  reduction  can  result  in  a  localized  mass  flow  reduction.  While  local  separation  is 
usually  present  on  the  blades  of  a  compressor,  the  extent  to  which  it  covers  a  blade  row 
only  becomes  relevant  if  there  is  a  corresponding  reduction  in  flow  or  a  coalescence  of 
many  adjacent  localized  separation  zones  into  a  stall  cell  which  leads  to  full  compressor 
stall. 
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A. 


TYPES  OF  STALL 


The  stall  phenomenon  discussed  in  this  document  can  be  separated  into  two  major 
categories:  Rotating  Stall  and  Compressor  Surge.  While  both  major  categories  can  trace 
their  roots  to  deviation  in  flow  properties  from  design  (Angle  of  Attack,  Density, 
Temperature,  etc.)  their  manifestation  and  effect  on  the  operation  of  a  compressor  can  be 
markedly  different. 

1.  Rotating  Stall 

Rotating  Stall  is  caused  by  the  coalescence  of  a  few  adjacent  locally  stalled  blade 
passages  into  a  cell(s)  that  rotates  in  the  same  direction  as  the  rotor  but  at  a  slower  speed, 
Figure  3.  This  rotation  is  caused  because  of  the  reduction  in  flow  rate  associated  with  a 
stall  cell(s).  As  blades  immediately  adjacent  to  a  stall  cell  begin  to  increase  their  flow 
rate  to  accommodate  the  cells  local  reduction,  they  experience  a  change  in  angle  of  attack 
that  directly  leads  to  separation  and  stall.  In  this  way,  a  stall  cell  can  “walk”  around  the 
annulus  of  a  compressor  rotor. 


Figure  3.  Rotating  Progressive  Stall  Cell 
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Rotating  Stall  can  manifest  itself  as  a  stalled  region  extending  the  full  span  of  the 
blade  (abrupt  stall)  or  only  a  partial  span  of  the  blades  (progressive  stall).  In  both  cases, 
the  stall  cell(s)  generally  rotate  between  50-70%  of  the  blade  speed  and  are  characterized 
by  excess  noise,  vibration,  and  a  reduction  in  performance.  The  primary  difference  is  in 
the  magnitude  of  the  disturbance  to  the  system.  In  the  case  of  progressive  stall 
performance,  degradation  is  usually  minor  and  the  only  external  indication  of  stall  is  an 
increase  in  noise  levels  and  vibration.  In  abrupt  stall,  all  of  these  symptoms  are 
magnified  with  the  potential  of  the  stall  cell(s)  extending  through  multiple  compressor 
stages  and  a  significant  reduction  or  total  loss  of  pressure  rise  across  the  compressor 
resulting  in  surge. 

2.  Compressor  Surge 

While  Rotating  Stall  is  a  phenomenon  associated  with  the  interaction  of  adjacent 
blades  on  a  rotor,  Compressor  Surge  is  due  to  an  unstable  fluid  coupling  of  the  rotor  and 
stator.  This  is  often  characterized  by  a  violent  oscillation  in  rotor  speed. 

The  unstable  interaction  between  the  rotor  and  stator  is  due  to  variation  in  fluid 
flow  properties  from  design  parameters.  Similar  to  rotating  stall  a  cell(s)  is  formed; 
however,  in  the  case  of  the  Sanger  Stage,  it  is  thought  to  that  the  stator  blades  are  first 
affected.  The  stalled  stator  chokes  the  flow  locally  and  reduces  mass  flow  through  the 
rotor  blades,  stalling  and  chocking  the  flow  through  the  rotor,  and  reducing  mass  flow  to 
the  stator.  This  in-tum  allows  the  stator  to  recover  and  operate  normally.  If  the 
underlying  cause  of  the  stall  is  not  corrected  then  this  stall-recovery  cycle  will  repeat, 
Figure  4. 
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Figure  4.  Compressor  Surge  Cycle 


B.  IMPACT  OF  COMPRESSOR  STALL  OR  SURGE  ON  GAS  TURBINES 
ENGINES 

Gas  turbine  engines  are  high-precision  machines  designed  to  run  at  steady-state 
conditions  around  predetermined  operating  points.  Since  blade  stall  manifests  itself  as 
highly  disturbed  flow  (generating  rapid  changes  in  pressure,  velocity,  and  fluid  density)  it 
can  impact  the  performance  of  a  turbine  in  numerous  ways  ranging  from  minor 
degradation  in  thrust  and  power  to  a  complete  mechanical  engine  failure. 

1.  Performance  Reduction 

Compressor  stall  at  its  most  basic  level  represents  a  reduction  in  intake  airflow 
through  a  gas  turbine.  This  reduction  in  flow  through  the  engine,  depending  on  its 
severity,  can  result  in  a  measureable  performance  reduction  manifesting  as  reduction,  or 
loss,  in  thrust  and/or  a  reduction  in  fuel  economy. 

If  the  stall  cell  grows  into  a  surge  cycle  it  may  effectively  stop  the  flow  of  high- 
pressure  air  to  the  burner  section.  The  resulting  inability  of  the  compressor  blades  to 
impart  mechanical  work  to  the  fluid  can  allow  the  high-pressure  fluid  in  the  burner 
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section  to  reverse  direction  and  exit  the  intake  section  of  the  engine  or  allow  excess 
unbumed  fuel  to  be  expelled  from  the  exhaust  section.  In  both  cases,  when  fuel  rich  air  is 
allowed  to  exit  the  turbine,  it  may  result  in  secondary  burning  of  the  fuel  or  an  explosion 
when  it  is  allowed  to  mix  with  oxygen-rich  air  (as  shown  in  Figure  1). 

2.  Mechanical  Damage 

The  cyclical  nature  of  the  stall  cell,  as  it  “walks”  around  the  compressor  annulus 
or  progressively  stalls  and  recovers,  results  in  fatigue  loading  well  in  excess  of  engine 
design  specifications.  The  asymmetrical  nature  of  the  compressor  loading  under  stall 
also  imparts  unsteady  loads  to  the  bearings  and  support  structures  designed  primarily  for 
thrust  loads.  Both  effects  work  to  decrease  engine  life  and  increase  lifecycle  maintenance 
costs. 


3.  Recovery 

Once  a  compressor  has  entered  surge  it  will  not  recover  unless  there  is  a  change  in 
airflow  properties  (generally  produced  by  a  change  in  throttle).  The  change  in  throttle 
required  to  effect  recovery  is  called  compressor  hysteresis.  Figure  5  shows  compressor 
surge  taken  at  a  constant  throttle  setting.  The  level  of  throttle  reduction,  hysteresis, 
required  to  return  to  steady  state  conditions  is  shown.  Once  corrective  action  is  taken,  the 
compressor  returns  to  a  state  that  is  not  the  same  as  the  point  where  the  compressor 
entered  the  surge  cycle. 

The  magnitude  of  the  hysteresis  is  highly  dependent  on  compressor  design 
increasing  as  the  pressure  ratio  across  stages  increases;  however,  large  pressure  ratios  are 
extremely  desirable  because  they  represent  an  increase  in  engine  efficiency  and 
performance  as  well  as  a  reduction  in  weight  and  size  for  a  given  compressor.  These 
competing  design  requirements  can  manifest  in  lower  pressure  ratios  or  in  higher  control 
system  costs  to  help  reduce  the  amount  of  time  required  to  restart  a  stalled  engine  in 
flight. 
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Hysteresis 


Figure  5. 


Surge  Cycle  Showing  Reverse  Flow  and  Hysteresis,  After  [12] 
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III.  OVERVIEW  OF  PREVIOUS  RESEARCH 


Initial  research  conducted  on  the  transonic  compressor  rig  focused  on  establishing 
the  stall  margins  at  various  speed  (Table  2)  and  air  only  inlet  flow  conditions.  Steady 
State  and  Transient  analysis  of  the  rig  in  its  rotor-only  configuration  led  to  the 
verification  of  these  parameters  as  well  as  identification  of  low  frequency  instabilities 
prior  to  and  during  stall.  Later  on  additional  rotor-only  analysis,  to  identify  reductions  in 
stall  margin,  with  high  and  low  pressure  steam  ingestion  was  also  conducted. 

Research  using  the  full  stage,  rotor  and  stator,  was  conducted  to  identify  changes 
in  performance  characteristics  and  stall  margin  when  compared  to  the  rotor  only 
configuration  in  both  throttle  and  steam  induced  stall. 


A.  ROTOR  ONLY  ANALYSIS 

Rodgers  reestablished  measurements  and  testing  procedures  at  60%,  70%,  and 
80%  for  the  Sanger  developed  transonic  compressor  rig  at  the  TPL  [2],  This  research 
was  expanded  by  Villescas  [3],  Brunner  [4],  and  Payne  [5]  to  include  full  pressure 
measurements  at  60%,  70%,  80%,  90%,  95%,  and  100%  of  rotor  speed,  Figure  6.  Steady 
and  unsteady  analysis  of  pressure  measurements  and  identification  of  stall  margins  prior 
to  steam  introduction  were  presented. 


U 

i 


-  tett: 

4  4  ■*  bMW 

-  tew  t.iN 

m 

j*<b 

1 

hi 

hhJ 

*  k- 

1 

-  \ 

#0.  *1 

- 4 

k 

W 

- - 3 

1 

l 

k 

*  * 

K 

T 

A 

mi 

ki 

<  <4  *  t  1»  T  T*  «  St 

*-— *'"lfin**i 


Figure  6.  Established  Stall  Margin  and  Pressure  Mapping,  From  [5] 
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Villescas  [3]  further  investigated  rotor  diffusion  factors  at  choke,  peak  efficiency, 
and  stall  as  a  means  to  predict  compressor  stall.  Direct  measurement  and  analysis 
showed  a  clear  correlation  between  diffusion  factor  and  the  onset  of  compressor  stall, 
Figure  7. 


Figure  7.  Diffusion  Factor  at  Near  Stall,  From  [3] 

Flow  field  properties  were  also  characterized  to  facilitate  future  Computation 
Fluid  Dynamics  (CFD)  modeling  with  measurements  of  inlet  and  outlet  Mach  number 
(Figure  8)  and  Flow  Angle  (Figure  9)  as  a  function  of  the  hub-tip  radius  ratio,  Villescas 
[3]  and  Brunner  [4], 
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Figure  9.  Inlet  Flow  Angle  Measurements  at  100%  Rotor  Speed,  From  [3] 

The  Pitch  Angle  at  the  compressor  inlet  was  measured  at  the  peak  efficiency 
points  at  70%,  80%,  and  90%  rotor  speeds,  Figure  10.  This  data  identified  inlet  flow 
swirl  at  the  inlet  of  the  compressor  further  characterizing  the  flow  for  CFD  analysis. 


Figure  10.  Magnitude  of  Flow  Pitch  Toward  the  Center  Line,  From  [4] 


A  CFD  model  of  the  flow  inlet  ducting  and  bell-mouth  were  generated  using  the 
Ansys  CFX-5  software  package,  Figure  11. 
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Figure  1 1 .  Computation  Fluid  Dynamics  Model  of  Compressor  Inlet,  From  [4] 


Computational  results  were  then  compared  to  the  experimental  data  as  a  means  of 
validating  the  model  for  future  redesign  and  analysis,  Figure  12.  A  strong  correlation 
between  the  numerical  and  experimental  results  allowed  for  future  modification  of  the 
experimental  apparatus  prior  to  steam  induction  testing. 


Figure  12.  Inlet  Mach  Number  Experimental  Data  and  CFD  Comparison,  From  [4] 


A  full  analysis  of  efficiency  and  pressure  ratio  under  varying  rotor  speeds  was 
presented  by  Payne  [5],  Using  Fast  Fourier  Transform  (FFT)  Techniques  to  determine 
the  dominant  frequencies  and  their  magnitudes  stall  precursors  were  identified  in  the  pre¬ 
stall  condition,  Figure  13. 
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Figure  13.  Power  Spectrum  at  90%  Speed  Near  Stall,  From  [5] 

Evaluating  the  power  spectrum  data  at  fixed  intervals  in  time  verified  that  the 
precursor  modal  instabilities  were  stable  and  present  throughout  the  sampling  period, 
Figure  14.  This  plot  also  displayed  a  clear  boundary  between  the  flow  properties 
immediately  preceding  and  following  stall.  The  location  of  the  stall  cell  at  roughly  60% 
of  the  blade  passage  frequency  with  no  indication  of  instability  at  any  frequency 
harmonics  indicates  a  single  stall  cell  rotating  around  the  compressor  face. 
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Figure  14.  Power  Contour  Plot  of  Kulite  Data  at  70%,  From  [5] 
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Using  hotfilm  data  at  varying  speeds  and  mass  flow  rates,  an  average  turbulence 
intensity  was  calculated  and  plotted.  While  no  clear  correlation  has  found  between  the 
rotor  speeds,  mass  flow  rate,  or  turbulence  intensity,  an  average  intensity  of  2.26%  was 
determined,  Figure  15.  This  further  characterization  assisted  future  CFD  modeling. 


Figure  15.  Turbulence  Intensity  vs.  Mass  Flow  Rate,  From  [5] 


A  CFD  model  of  the  compressor  rotor  was  then  generated  using  Ansys  CFX-5, 
Figure  16.  To  reduce  the  computational  requirements,  a  single  blade  passage  was 
modeled  with  periodic  boundary  conditions  to  simulate  a  full  rotor. 
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Computation  Fluid  Dynamics  Model  of  Compressor  Rotor,  From  [6] 
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Figure  16. 


This  model  was  run  at  varying  speeds  and  inlet  conditions,  and  the  results  were 
compared  to  the  experimentally  derived  results,  Figure  17.  The  model  successfully 
approximated  the  rotor  performance  characteristics  with  an  underestimation  of  the 
pressure  ratio  and  overestimation  of  the  efficiency.  This  minor  discrepancy  indicated 
that  a  more  complete  understanding  of  inlet  conditions  was  required  prior  to  the 
generation  of  a  new  CFD  model. 


Figure  17.  CFD  to  Experimental  Results  Comparison  of  Pressure  Ratio,  From  [6] 


Work  accomplished  by  Koessler  adjusted  the  compressor  stall  margin  using  low- 
and  high-pressure  steam  ingestion,  Figure  18.  An  adjustment  of  the  stall  margin  of  up  to 
30%  was  identified. 
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B.  FULL  STAGE  ANALYSIS 


Hurley  established  the  performance  map  for  the  full  compressor  stage  at  70%, 
80%,  90%,  95%,  and  100%  of  rated  speed  using  throttle  and  steam  stall  procedures, 
Figure  19.  He  identified  a  reduction  in  stall  margin  for  the  throttled  stage  when 
compared  to  a  steam  induced  stall  from  3-5%. 


Figure  19.  Compressor  Stage  Performance  Map,  From  [1 1] 

A  comparison  of  the  stall  margins  of  the  rotor-only  and  stage  configurations 
identified  a  significant  reduction  in  margin  of  10-14%  under  steam  and  throttle  induced 
stall,  Figure  20.  This  overall  reduction  in  performance  across  all  rotor  speeds  and  throttle 
settings  indicated  that  the  stator  was  entering  stall  prior  to  the  rotor. 
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Figure  20.  Compressor  Stage  and  Rotor-Only  Performance  Comparison,  From  [1 1] 

An  interpolation  of  this  pressure  data  resulted  in  pressure  contour  plots  that 
clearly  showed  oblique  shocks  at  the  leading  edge  of  the  rotor  blades,  Figure  21.  This 
visual  representation  of  the  flow  field  allows  a  refinement  and  validation  of  past  CFD 
models. 


Figure  21 .  Rotor-Only  Pressure  Contour  Mapping,  From  [11] 

A  continuation  of  this  research  will  be  presented  in  the  following  chapters. 
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IV.  EXPERIMENTAL  APPARATUS 


A.  TRANSONIC  COMPRESSOR  RIG 

Unlike  previous  work  that  focused  only  on  one  test  rig  configuration,  this  research 
employed  data  collected  from  test  runs  in  the  rotor-only  and  full  stage  configuration  of 
the  Sanger  test  stage  [1].  This  experimental  apparatus  has  been  described  in  detail  in 
previous  research  and  is  only  briefly  summarized  below. 


Figure  22.  Sanger  Test  Stage  in  Rotor-only  and  Full  Stage  Configuration 

The  full  stage  employed  at  the  NPS  Turbo  Propulsion  Laboratory  consisted  of  a 
22-blade  rotor  and  a  27-blade  stator,  Figure  22.  Both  were  manufactured  from  T7075-T6 
Aluminum  alloy  because  of  its  lightweight  and  durability.  The  later  was  particularly 
important  because  of  the  strong  vibration  and  high  possibility  of  mechanical  damage 
during  the  stall  events  that  culminated  each  speed  run.  A  detailed  description  of  the 
critical  design  parameters  is  shown  in  Table  1. 
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Rotor  Pressure  Ratio 

1.61 

Stage  Pressure  Ratio 

1.56 

Tip  Speed 

396.2  m/sec 

Design  Speed 

27,085  rpm 

Design  Mass  Flow 

7.75  kg/sec 

Specific  Mass  Flow 

170.9  kg/sec-m2 

Specific  Head  Rise 

0.246 

Tip  Inlet  Relative  Mach  Number 

1.6 

Aspect  Ratio 

1.2 

Hub/Tip  Radius  Ratio 

0.51 

Number  of  Rotor  Blades 

22 

Number  of  Stator  Blades 

27 

Tip  Solidity  (Rotor) 

1.3 

Tip  Solidity  (Stator) 

1.0 

Outside  Diameter 

0.2794  m  (11”) 

Rotor  Diffusion  Factor  (Tip) 

0.40 

Rotor  Diffusion  Factor  (Hub) 

0.47 

Stator  Diffusion  Factor  (Tip) 

0.52 

Stator  Diffusion  Factor  (Hub) 

0.58 

Table  1.  Sanger  Test  Stage  Design  Parameters 


Rotor  speed  was  controlled  through  two  single-stage  power  turbines,  on  a 
common  shaft,  that  received  air  from  a  12-stage  Allis-Chalmers  Axial  compressor  that 
took  air  directly  from  the  atmosphere  and  delivered  it  directly  to  the  TCR  at  the  required 
flow  rate.  The  compressed  air  used  to  drive  the  axial  turbine  was  throttled  to  the  flow 
rate  needed  to  maintain  the  desired  rotor  speed  and  was  exhausted  back  to  atmospheric 
radially.  Flow  through  the  compressor  stage  was  controlled  by  an  upstream  throttle  plate 
that  receives  electrical  command  from  the  test  control  station. 
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Figure  23.  Transonic  Test  Rig  Schematic,  After  [11] 


As  discussed  by  Hurley  the  test  rig  differs  from  previous  research  due  to  the 
removal  of  a  honeycomb  flow  straightener  located  behind  the  fan  and  the  introduction  of 
ceramic  high  temperature  ball  bearings  [11]. 

B.  DATA  ACQUISITION  SYSTEM 

Data  collection  was  accomplished  using  a  high-speed  and  low-speed  system  for 
recording  pressure  and  temperature  measurements  across  the  rotor  as  well  as  rotor  speed 
measurements. 

1.  Slow-Speed  Data  Collection 

Steady-state  measurements  were  taken  using  pressure  probes  positioned  upstream 
of  the  rotor  and  at  the  same  axial  position  of  the  Kulite  probes  so  stagnation  pressure 
measurements  could  be  taken  for  each  variation  of  rotor  speed  and  throttle  position. 
Additional  measurements  were  also  made  via  1 1  stagnation  temperature  probes 
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positioned  downstream  of  the  rotor.  Using  these  probes  isentropic  efficiency  and  mass 
averaged  temperature/pressure  ratios  were  calculated  and  use  to  determine  the 
compressor  performance  maps. 


2.  High-Speed  Data  Collection 
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Figure  24.  Schematic  of  High  Speed  Data  Capture  System,  After  [11] 

As  shown  in  Figure  23  and  Figure  24,  the  high-speed  system  consists  of  ten 
Kulite  pressure  probes  in  the  casewall  across  the  rotor  blades,  one  Kulite  impact  probe 
upstream,  and  two  hot  film  probes.  The  Kulite  probes  are  capable  of  sampling  at  300 
kHz  and  the  hotwire  probes  can  respond  to  changes  in  velocity  at  400  Hz.  They  are 
sampled  and  stored  via  a  VXI  Mainframe  at  its  maximum  sampling  frequency  of  196.6 
kHz  to  reduce  the  aliasing  of  high-frequency  signals  into  the  lower  frequency  regions  of 
interest.  Data  is  stored  in  a  binary  format  internally  on  the  VXI  Mainframe  but  can  be 
exported  to  PC  in  a  delimited  text  file  for  analysis.  This  file  format  is  used  to  ensure 
portability  of  data  regardless  of  system  type  or  post  processing  software. 
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As  shown  in  Figure  25  the  Kulite  probes  are  positioned  on  the  casewall  so  that 
static  pressure  measurements  can  be  taken  across  the  chord  of  the  rotor  blades  starting  at 
a  position  well  upstream  (Probes  1 ,  7  and  8)  and  ending  at  a  position  well  downstream 
(Probe  6).  The  upstream  probes  are  positioned  120-degrees  around  the  casewall  at  the 
same  axial  location  so  that  the  speed  of  any  stall  cells  can  more  accurately  be  determined. 
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V.  EXPERIMENTAL  PROCEDURE 


As  in  previous  research  conducted  on  the  Transonic  Compressor  Rig  (TCR) 
steady  state  reading  were  recorded  using  the  high  and  low  speed  data  so  that  a 
compressor  performance  map  could  be  generated.  This  was  accomplished  by  running  the 
compressor  at  a  constant  RPM  and  varying  throttle  position  so  that  the  progression  of 
pressure  ratio  and  efficiency  could  be  recorded  and  plotted,  Figure  26. 
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Figure  26.  Full  Stage  Compressor  Performance  Maps 


To  obtain  the  data  from  the  low-  and  high-speed  data  acquisition  system,  the 
compressor  rig  was  operated  in  steady  state  condition  at  one  of  the  speeds  identified  in 
Table  2.  This  speed  was  then  corrected  for  local  variation  in  temperature  and  pressure  so 
that  the  inlet  Mach  number  for  any  given  run  would  be  the  same  regardless  of  daily 
atmospheric  variation. 


%  Rotor  Speed 

RPM 

100% 

27,085 

95% 

25,730 

90% 

24,375 

80% 

21,670 

70% 

18,960 

Table  2.  Rotor  Speed  Profile 
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A  generic  test  procedure  consists  of  the  following: 

•  The  rotor  RPM  is  stabilized  around  the  corrected  RPM  for  the  desired  speed 
using  the  Allis-Chalmers  Controls. 

•  The  slow-speed  data  from  the  stagnation  pressure  and  temperature  probes  are 
recorded  via  a  PC  based  data  acquisition  system. 

•  The  high-speed  data  is  recorded  using  varying  back  pressures  on  the  Kulite 
Probes  (0,  5,  10  and  15psi).  This  is  done  so  that  the  data  can  be  accurately 
calibrated  and  analyzed.  For  steady  state  operation,  a  2-second  recording  is 
taken  using  the  high-speed  system;  for  a  stall  run,  a  60-second  recording  is 
taken  to  ensure  that  the  full  stall  event  is  captured. 

•  The  throttle  is  then  set  to  a  new  position  in  preparation  of  the  next  run. 
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VI.  DATA  ANALYSIS  DESCRIPTION 


The  underlying  MATLAB  code  used  in  this  research  represents  a  revision  of  the 
code  primarily  developed  by  A.  J.  Gannon  and  has  been  used  in  most  previous  theses  on 
the  TCR.  The  contribution  of  this  research  was  an  improved  implementation  of  the  data 
loading  and  calibration  routines  and  the  filter  construction  and  application  algorythm.  A 
simple-to-use  graphical  user  interface  (GUI)  was  also  developed  to  enable  rapid  data 
analysis. 

A.  USER  INTERFACE  FUNCTIONALITY 


Figure  27.  Top  Level  GUI  Functionality 


The  Data  Analyzer  GUI  was  used  to  accomplish  five  primary  tasks,  Figure  27: 

1)  Load  and  calibrate  raw  voltage  signals  from  the  Kulite  pressure 
transducers. 

2)  Perform  a  steady-state  analysis  of  the  pressure  signals,  in  the  time  and 
frequency  domains,  to  identify  precursor  frequencies  and  bulk  inlet  flow 
properties. 

3)  Generate  a  time-averaged  steady-state  pressure  contour  of  the  blade 
passages.  This  enables  the  identification  of  flow  field  structure  within  a 
typical  blade  passage. 
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4)  Generate  a  Transient  Pressure  Contour  plot  of  all  the  rotor  blades  over  a 
given  set  of  rotations  to  visualize  the  flow  field  as  it  changes  in  the 
passages  over  time. 

5)  Perform  a  Transient  Contour  Analysis  of  the  pressure  data  at  discrete 
locations  within  a  blade  passage.  Information  regarding  magnitude  and 
frequency  and  standard  deviation  of  the  locations  are  presented. 

B.  RAW  DATA  CALIBRATION 

1.  MATLAB  Code  Overview 

This  Data  Loading  and  Calibration  subroutine  used  the  m-files  briefly  described 
below  (Figure  28)  to  generate  the  binary  files  used  in  all  subsequent  analysis. 


Figure  28.  Data  Load  and  Calibrate  Workflow 


2.  Calibration  Description 

The  raw  Kulite  data,  recorded  using  varying  back  pressures,  was  processed  using 
the  loadkulitedata.m  script  so  that  it  could  be  accurately  converted  from  voltages  to 
pressure  ratios  with  regard  to  the  recorded  atmospheric  pressure  for  that  particular  run. 
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The  average  of  the  raw  voltage  for  each  Kulite,  the  static  pressures  recorded  from  probes 
at  the  same  axial  location  on  the  compressor  casing,  and  the  pressure  taken  from  the 
impact  Kulite  positioned  upstream  of  the  rotor  were  processed.  These  calibration  curves, 
uniquely  determined  for  each  Kulite,  were  then  applied  to  all  voltages  recorded  during  a 
given  series  of  test  runs  and  were  used  to  develop  all  the  experimental  results  discussed 
in  subsequent  chapters,  Figure  29. 


Figure  29.  Generic  Raw  Kulite  Voltage  Conversion  to  Pressure  Ratio  Data  Sets 

C.  STEADY  STATE  ANALYSIS  ROUTINE 
1.  MATLAB  Code  Overview 

This  Steady-State  Analysis  subroutine  used  the  m-files  briefly  described  below 
(Figure  30)  to  perform  the  analysis  and  generate  the  data  plots. 
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Figure  30.  Steady  State  Analysis  Workflow 


2.  Analysis  Description 

The  steady-state  analysis  was  conducted  by  transferring  the  pressure  data  into  the 
frequency  domain  so  that  frequency  regions  of  interest  could  be  located  and  data 
processing  methods  could  be  better  targeted,  Figure  31.  All  frequencies  were  normalized 
with  the  blade  passage  frequency  at  the  given  rotor  speed  so  that  the  dominant  blade 
passing  frequencies  could  be  more  easily  identified  at  their  integer  locations.  Primary 
areas  of  interest  were  frequencies  lying  below  and  just  above  the  blade  passing 

frequency  (l0° )  because  these  regions  were  known  locations  of  low-frequency  oscillation 

at  all  rotor  speeds.  However,  the  relative  magnitude  of  the  multiples  of  the  blade  passing 
required  that  a  filter  be  developed  so  that  the  relatively  small  peaks  within  the  regions  of 
interest  could  be  more  easily  identified. 
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Figure  3 1 .  Fast  Fourier  Transform  of  Pressure  Signal  at  1 00%  Rotor  Speed 

To  construct  the  filter  built-in  MATLAB  Fast  Fourier  Transform  and  filtering 
functions  were  employed.  A  Blackman  Window  was  selected  as  the  base  filter  because 
of  its  high  roll-off  and  low  ripple,  Figure  32. 


Figure  32.  Generic  Blackman  Window  Employed  in  Analysis 
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The  Blackman  window  was  sized  so  that  the  signal  would  be  attenuated  as 
quickly  as  possible  outside  of  the  regions  of  interest  without  adversely  impacting  the  raw 
signal  in  the  frequency  range  of  interest.  The  exact  size  was  determined  using  the 
sampling  frequency  (  /^ )  and  the  desired  roll-off. 
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A  roll-off  of  1 -Decade  was  ultimately  chosen  because  of  the  quick  attenuation  it 
provided  to  signals  immediately  outside  of  the  analysis  region  and  because  it  limited  the 
window  size  to  a  manageable  length.  The  filter  elements  were  added  together  to  generate 
arbitrary  band  pass  filters  allowing  an  analysis  across  any  frequency  range,  Figure  33. 
This  also  had  the  benefit  of  providing  a  reduction  of  high  frequency  sensor  noise  and  low 
frequency  sampling  error. 


Figure  33.  Signal  with  Band-Pass  Filtering  Applied 


This  frequency  data  was  then  used  to  generate  the  waterfall  plots  of  the 
oscillations  across  the  blade  passage  as  a  function  of  frequency  and  percent  of  chord. 
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The  data  is  then  transferred  back  into  the  time  domain  using  MATLAB’s  inverse  Fourier 
Transform  function  so  that  bulk  inlet  properties  can  be  calculated. 

D.  STEADY  STATE  CONTOUR  ROUTINE 
1.  MATLAB  Code  Overview 

This  Steady-State  Pressure  Contour  subroutine  used  the  m- files  briefly  described 
below  (Figure  34)  to  perform  the  analysis  and  generate  the  data  plots. 


Figure  34.  Steady  State  Pressure  Contour  Workflow 

2.  Analysis  Description 

The  Steady-State  Contour  program  used  the  unfdtered  Kulite  data  to  generate  a 
snapshot  of  a  time-averaged  typical  blade  passage.  Due  to  the  fact  that  pressure  data  only 
exists  at  the  angular  Kulite  positions  (The  vertical  lines  from  Figure  35)  the  axial  pressure 
data  had  to  be  interpolated  using  MATLAB.  This  required  that  all  Kulite  data  be  lined  up 
as  if  they  existed  along  a  single  angular  location.  The  data  was  then  interpolated  in  three 
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distinct  regions  to  reduce  the  interpolation  error  in  locations  of  known  high  gradient 
(Along  the  shock  (Rl),  horizontally  through  the  blade  passage  (R2),  and  along  the  blade 
(R3).)  These  regions  were  then  stitched  together  using  modified  cosine  functions  so  that 
the  transitions  were  seamless. 
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Figure  35.  Steady  State  Pressure  Contour  Plot  Showing  Interpolation  Regions 

E.  TRANSIENT  CONTOUR  ROUTINE 
1.  MATLAB  Code  Overview 

This  Transient  Pressure  Contour  subroutine  used  the  m-files  briefly  described 
below  (Figure  36)  to  perform  the  analysis  and  generate  the  data  plots. 
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Figure  36.  Transient  Pressure  Contour  Workflow 

2.  Analysis  Description 

This  program  performed  the  contouring  interpolation  in  the  same  manner  as  the 
Steady-State  Contour  sub-routine  with  the  exception  that  it  did  not  perform  the  time 
averaging.  It  used  the  high-pressure  point  at  the  leading  edge  of  the  blade  to  identify  the 
blade  tip  and  assumed  a  constant  pitch  length  to  correctly  interpolate  the  shock  and  blade 
angle.  In  addition  to  the  time  strip  plots  of  the  identified  rotations  the  program  also 
generated  a  video  file  by  using  individual  blade  passages  as  still  frames.  The  file  was 
then  saved  as  a  compressed  AVI  file  or  a  high-quality  uncompressed  AVI  file. 

F.  TRANSIENT  CONTOUR  ANALYSIS  ROUTINE 

1.  MATLAB  Code  Overview 

This  Transient  Contour  Analysis  subroutine  used  the  m-files  briefly  described 
below  (Figure  37)  to  perform  the  analysis  and  generate  the  data  plots. 
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Figure  37.  Transient  Contour  Analysis  Workflow 

2.  Analysis  Description 

The  underlying  code  implemented  in  this  analysis  routine  was  based  on  the  same 
methods  employed  in  the  transient  contour  generation.  This  method  differed  in  that  it 
used  the  tach  signal  and  pressure  signal  to  locate  the  blade  tip  and  pitch  length  in  the 
pressure  data.  Because  the  sampled  pressure  data  did  not  sample  the  same  position  in  the 
passage  at  each  blade  passing  the  data  was  interpolate  at  regular  positions  so  that  the 
frequency  analysis  and  statistical  analysis  were  being  used  properly.  The  pressure  at  each 
discrete  location  was  the  processed  using  MATLAB’s  Fourier  Transfer  and  Standard 
Deviation  functions  so  that  the  regions  of  highest  change  would  be  identified. 
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VII.  ANALYSIS  AND  RESULTS 


This  research  represents  an  analysis  of  the  Transonic  Compressor  Rig  at  70%, 
80%,  90%  and,  100%  of  rated  speeds  in  both  the  rotor-only  and  full  stage  configuration 
to  better  quantify  the  differences  between  both  operating  conditions.  For  the  steady-state 
operating  mode,  an  analysis  of  the  raw  pressure  data  was  conducted  in  the  frequency 
domain  with  the  intent  of  determining  dominant  frequencies  and  the  effect  that  the 
addition  of  a  stator  had  on  the  dynamic  of  the  system  when  compared  to  previous  rotor- 
only  research.  A  separate  analysis  was  conducted  in  the  time  domain  with  the  intent  of 
correlating  the  physical  structure  of  the  flow  field  (Shocks,  Boundary  Layers,  etc.)  with 
the  localized  the  frequency  domain  data.  A  transient  analysis  of  the  TCR  as  it  entered 
stall  was  then  conducted  using  many  of  the  same  tools  and  methods  developed  for  steady 
state  operation. 

A.  FREQUENCY  DOMAIN 

Using  the  filtered  Kulite  pressure  signals  described  in  Chapter  IV,  it  became 
possible  to  display  the  power-frequency  data  across  the  compressor  blade  so  that  the 
relative  values  and  frequencies  could  be  used  to  identify  frequency-spatial  regions  of 
interest  for  further  analysis  in  the  time  domain.  The  previously  described  frequency 
normalization  was  used  in  these  plots;  however,  the  Axial  Location  was  also  normalized 
against  the  blade  chord  so  that  “0”  denotes  the  leading  edge  and  “100”  denotes  the 
trailing  edge,  Figure  38. 
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Figure  38.  Full  Stage  Waterfall  Plot  of  Filtered  Frequency  Data  at  80%  Speed, 

Near  Stall 


Using  the  waterfall  plot  of  the  filtered  frequency  data  for  all  Kulite  probes  the 
frequency  location,  peak  magnitude,  and  extent  of  blade  coverage  can  be  determined  for 
any  given  speed-throttle  profile.  The  data  displayed  for  80%  speed  in  Figure  38  clearly 
shows  two  dominant  peaks  at  0.5  and  1.5  times  the  blade  pass  frequency  with  the  higher 
frequencies  dominating  the  blade  tip  and  lower  frequencies  more  prominent  at  the  trailing 
edge.  This  hints  at  the  presence  of  multiple  instabilities  rotating  around  the  front  of  the 
rotor  and  a  single  instability  in  front  of  the  stator. 

Waterfall  plots  were  generated  for  all  speed-throttle  settings  so  that  the 
progression  of  instability  growth  as  the  rotor  approached  stall  could  be  clearly  analyzed. 
Figure  39  shows  this  progression  for  the  TCR  at  100%  of  rotor  speed  in  the  rotor-only 
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configuration.  As  shown  in  Block  B,  strong  peaks  appeared  at  0.5  times  the  blade 
passage  frequency  as  the  rotor  approached  stall  and  suddenly  reduced  in  power  in  the  test 
run  just  preceding  the  onset  of  stall. 


Figure  39.  Rotor  Only  Waterfall  Plots  of  Filtered  Frequency  Data  at  100% 

Rotor  Speed 


Another  important  feature  of  these  plots  is  the  extent  of  the  blade  that  was 
experiencing  these  high-pressure  variations.  In  most  cases  there  was  a  substantial 
“wake”  trailing  from  the  peaks  experienced  at  the  leading  edge  of  the  blade;  however,  in 
all  the  tests  below  95%  speed  there  was  also  a  substantial  transport  of  information 
forward  of  the  blade,  Figure  40. 
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100%  Rotor  Speed  95%  Rotor  Speed  90%  Rotor  Speed 

Figure  40.  Waterfall  Plots  of  Filtered  Frequency  Data  Showed  Forward  Transport 

Figure  41  shows  the  same  analysis  for  the  95%  rotor  speed  run  in  the  full  stage 
configuration.  There  does  appear  to  be  a  series  of  peaks  whose  magnitude  and  location  is 
stable  in  the  vacinity  of  0.1  and  0.9  times  the  blade  pass  frequency;  however,  these 
locations  do  not  correspond  to  known  pre -modal  instabilities  and  are  usually  only  found 
at  the  location  of  a  single  Kulite  probe  hinting  at  an  unstable  shock  boundary  interaction 
being  registered. 
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Figure  41.  Full  Stage  Waterfall  Plots  of  Filtered  Frequency  Data  at  95%  Rotor  Speed 


As  a  means  to  quantify  the  total  power  present  in  the  rotor  at  a  given  speed- 
throttle  profile  the  RMS  power  was  calculated  in  the  frequency 

range {0. 1  < /  <  0.9, 1.1  </  <  1.9}  by  summing  the  peak  values  along  a  given  Kulite  and 


then  averaging  the  power  of  all  the  Kulites. 
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This  data  was  then  tabulated  and  plotted  for  the  full  stage  and  rotor  only 
configuration,  Figure  42.  As  shown  the  RMS  power  in  the  rotor  only  configuration 
increased  as  the  rotor  speed  increased  and  as  the  throttle  was  closed.  This  behavior  was 
not  present  in  the  stage  where  power  levels  remained  constant  even  as  the  throttle  began 
to  close. 
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Figure  42.  RMS  Power  of  Rotor  and  Stage  in  Low  Frequency  Band 


Further  investigating  the  relationship  between  the  RMS  power  and  the  onset  of 
stall  the  RMS  power  throughout  the  frequency  band  {0.1<  /  <44}  was  calculated  to 
include  the  blade  pass  frequency  peaks,  Figure  43. 
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Figure  43.  RMS  Power  of  Rotor  and  Stage  in  Full  Frequency  Spectrum 


The  total  frequency  RMS  plot  in  the  rotor-only  configuration  was  very  similar  in 
shape  when  compared  to  the  low-frequency  plots  only  differing  in  the  magnitude  of  the 
power  present.  This  is  contrasted  by  the  behavior  of  the  full  stage  which  now  began  to 
displayed  the  same  shape  and  magnitude  as  the  rotor  only  case.  This  prompted  the 
development  of  a  non-dimensional  form  of  the  RMS  power  so  that  the  differences  could 
be  more  easily  explained,  Figure  44. 
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Figure  44.  Non-Dimensional  RMS  Power  for  Rotor  and  Stage 


While  no  clear  correlation  was  found  between  the  RMS  power  at  any  frequency 
level  and  the  onset  of  compressor  stall  an  analysis  of  the  RMS  power  distribution 
revealed  clear  differences  between  the  rotor-only  and  stage  configurations.  As  shown  in 
Figure  43  the  total  power  of  the  system  was  generally  the  same  at  all  speed-throttle 
profiles;  however,  Figure  44  revealed  that  the  sonic  and  transonic  rotor  speeds  (80%  to 
100%)  displayed  a  downward  trend  in  the  stage  configuration  and  a  flat  trend  in  the  rotor 
only  configuration. 
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Table  3  represents  a  tabulation  of  the  RMS  powers  at  the  pre-stall  throttle  position 
for  each  speed.  As  shown  the  total  RMS  power  in  the  rotor  and  stator  were  roughly  equal 
while  the  greatest  differences  in  behavior  were  present  in  the  low  frequency  bands.  This 
hints  at  a  rotor-stator  coupling  that  had  the  affect  of  moving  the  power  to  the  higher 
frequency  regions  suggested  that  the  inclusion  of  the  stator  had  stiffened  the  system  and 
raised  its  natural  frequency.  In  addition,  the  system  was  being  damped  as  indicated  by 
the  reduction  in  magnitude. 
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Low  Freq 
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70% 

1.13E-02 
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2.21E-03 
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3.23E-02 
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8% 
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24% 

90% 

6.81E-02 

5.32E-05 

7.81E-04 

1.00E-01 

3.32E-05 

3.32E-04 

19% 

23% 

40% 

95% 

1.01E-01 

9.56E-05 

9.47E-04 

1.10E-01 

4.34E-05 

3.95E-04 

4% 

38% 

41% 

100% 

9.42E-02 

1.63E-04 

1.73E-03 

1.03E-01 

6.57E-05 

6.38E-04 

4% 

43% 

46% 

Table  3.  RMS  Power  Comparison 


B.  TIME  DOMAIN  ANALYSIS 

The  filtered  signal  was  then  transferred  back  into  the  time  domain  through  the  use 
of  MATLAB’s  inverse  FFT  algorithm  so  that  bulk  flow  Mach  Numbers  and  Inlet  Angles 
could  be  calculated  for  the  compressor  inlet.  This  allowed  a  further  characterization  of 
the  inlet  flow  and  identified  key  differences  between  the  rotor-only  and  stage 
configurations. 

The  inlet  Mach  number  for  the  bulk  flow  was  calculated  using  the  ratio  of 
stagnation  to  static  pressure  (Pr)  and  the  ratio  of  specific  heats  (7).  A  constant  specific 

heat  was  assumed  for  the  flow;  this  was  reasonable  because  of  the  relatively  small 
changes  in  temperature  and  pressure  in  the  inlet  section  of  the  compressor. 
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Figure  45.  Raw  Bulk  Mach  Number  Data  at  100%  Rotor  Speed 

The  average  inlet  Mach  number  was  calculated  for  all  speed-throttle  settings  and 
plotted  versus  flow  rate,  Figure  46.  Both  plots  display  similar  behavior  in  the  transonic 
ranges  (90%  -  100%)  but  show  drops  in  Mach  number  in  the  subsonic  and  sonic  rotor 
speeds  (70%  and  80%). 
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Figure  46.  Inlet  Mach  of  Bulk  Flow  in  Rotor-only  and  Full  Stage 


The  peak-to-peak  variation  for  a  given  test  run  was  calculated  using  the 
instantaneous  maximum  and  minimum  Mach  number  values.  This  data  was  tabulated 
and  plotted  versus  mass  flow,  Figure  47.  Both  configurations  displayed  similar  behavior 
for  all  speed-throttle  settings  with  increased  in  the  peak-to-peak  magnitude  as  the  rotor 
approached  stall.  In  all  cases,  with  the  exception  of  the  80%  run,  this  variation  was  as 
high  as  50%  of  the  mean  Mach  number. 
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Figure  47.  Peak-to-Peak  Mach  Variation  of  Bulk  Flow  in  Rotor-only  and  Full  Stage 


As  shown  in  Table  4,  Mach  numbers  present  in  the  rotor-only  configuration  in  the 
pre-stall  speed-throttle  setting  were,  on  average,  15%  higher  than  those  of  the  full  stage. 
Also  shown  is  a  mean  10%  increase  in  peak-to-peak  variation  of  the  bulk  flow  for  the  full 
stage  in  these  same  pre-stall  conditions  at  all  rotor  speeds.  There  were  inconsistencies  in 
average  Mach  number  data  in  the  full- stage  (revealed  by  the  reduction  in  Mach  from  70% 
to  80%).  These  inconsistencies  could  not  be  readily  explained  using  this  analysis. 
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%  Difference 

Speed 

Average 

Variation 

Average 

Variation 

Average 

Variation 

70% 

0.22 

0.08 

0.16 

0.08 

15.8% 

0% 

80% 

0.25 

0.10 

0.11 

0.17 

38.9% 

25.9% 

90% 

0.31 

0.09 

0.25 

0.11 

10.7% 

10.0% 

95% 

0.35 

0.08 

0.31 

0.07 

6.1% 

6.7% 

100% 

0.37 

0.10 

0.30 

0.11 

10.4% 

4.8% 

Table  4.  Inlet  Bulk  Mach  Number  Comparison 


A  similar  analysis  was  also  conducted  by  calculating  the  bulk  flow  inlet  angle 
relative  to  the  rotor  blades  so  that  the  interaction  between  the  changes  in  rotor  speed  and 
inlet  Mach  number  could  be  better  determined.  The  Axial  Fluid  Velocity  (Cz)  is 

calculated  using  the  Mach  number  (M)  from  equation  3,  gas  constants (7, R),  and  static 
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Figure  48.  Raw  Bulk  Inlet  Angle  at  100%  Rotor  Speed 


Average  inlet  angle  and  peak-to-peak  variation  were  calculated  in  the  same  way 
that  the  average  and  peak-to-peak  Mach  was  calculated.  This  data  was  then  plotted 
versus  the  mass  flow  rate  for  each  throttle-speed  setting,  Figure  49  and  Figure  50. 
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Figure  49.  Mean  Relative  Inlet  Angle  of  Bulk  Flow  in  Rotor-only  and  Full  Stage 


The  progression  of  average  relative  inlet  angle  with  respect  to  the  rotor  blades 
was  very  nearly  identical  in  both  rotor-only  and  full  stage  configurations.  The  only 
difference  was  that  the  rotor  only  configuration  was  able  to  operate  at  lower  flow  rates 
prior  to  stall  which  extended  the  curves  further  to  the  left. 
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Figure  50.  Peak-to-Peak  Variation  of  Relative  Inlet  Angle  in  Rotor-only 

and  Full  Stage 

As  shown  in  Figure  50  curves  for  the  peak-to-peak  magnitude  for  the  two 
operating  modes  were  different  in  shape  and  magnitude.  This  was  contrasted  by  inlet 
Mach  number  which  behaved  similarly  in  both  test  rig  configurations.  The  shapes  of 
these  curves  were  similar  to  the  RMS  power  plots,  Figure  42  through  Figure  44,  where 
rotor  only  operations  revealed  a  strengthening  of  power  measurements  as  the  flow  rate 
was  throttled  and  full  stage  operations  were  very  nearly  flat. 
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An  analysis  of  the  relative  inlet  angles  was  then  conducted  to  determine  the 
instantaneous  maximum  relative  inlet  angle  because  this  angle  was  related  to  the  loading 
of  an  individual  blade,  growth  and  separation  of  the  boundary  layer,  and  the  inception  of 
stall  cells.  This  data  was  plotted  versus  mass  flow  for  all  throttle-speed  settings  and  is 
displayed  in  Figure  5 1 . 
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As  shown  in  Table  5,  the  analysis  of  bulk  relative  inlet  angle  revealed  that  the  test 
rig,  in  the  rotor-only  and  full  stage  configurations,  stalled  at  approximately  the  same 
maximum  and  mean  angles  for  all  speeds  hinting  at  a  stall  event  being  initiated  by  the 
rotor  vice  the  stator  which  previous  research  had  theorized  that  the  stator  could  stall  first 
at  the  hub.  Also  shown  is  a  generally  more  chaotic  flow  for  the  rotor-only  configuration 
in  a  manner  similar  to  what  has  been  previously  identified  in  the  RMS  power  analysis. 


Rotor 

Speed 

Rotor  Only 

Full  Stage 

%  Difference 

Mean 

Maximum 

Variation 

Mean 

Maximum 

Variation 

Mean 

Maximum 

Variation 

70% 

75.3 

77.1 

3.68 

75.7 

76.9 

1.91 

0.26% 

0.13% 

31.7% 

80% 

75.1 

76.9 

3.17 

77.8 

79.9 

3.84 

1.77% 

1.9% 

9.56% 

90% 

72.5 

74.6 

3.47 

73.5 

74.9 

2.19 

0.68% 

0.20% 

22.6% 

95% 

71.1 

72.4 

2.70 

71.4 

72.2 

1.65 

0.21% 

0.14% 

24.1% 

100% 

70.2 

72.1 

3.79 

71.4 

72.6 

2.28 

0.85% 

0.35% 

24.9% 

Table  5.  Relative  Inlet  Angle  Comparison 


C.  PRESSURE  CONTOUR  ANALYSIS 

Pressure  contour  plots  of  the  blade  passage  were  generated  using  the  methodology 
previously  described  in  Chapter  IV.  This  enabled  the  visualization  of  the  flow  field 
around  the  blades  at  steady-state  operations  by  averaging  the  pressure  data  at  each  blade 
pitch  into  a  single  frame. 
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(a)  Throttle  Nearly  Open  /  Choked  Flow  (b)  Approaching  Stall 

Figure  52.  Rotor  Only  Pressure  Contour  Plot  at  100%  Rotor  Speed 


An  Analysis  of  these  pressure  contour  plots  revealed  several  structures  that  could 
be  indentified  and  used  to  explain  the  frequency  analysis  results  previously  shown. 
Figure  52a  and  Figure  52b  represent  flow  fields  on  opposite  ends  of  the  progression  from 
fully  open  throttle  to  the  pre-stall  throttle  setting.  These  two  were  selected  for  this 
sample  analysis  because  they  best  represented  the  types  of  structures  found  at  all  speed- 
throttle  profiles. 

As  shown  in  Figure  52a,  the  open  throttle  mode  of  operation  displayed  attached 
and  oblique  shocks.  The  oblique  shock  reflected  off  the  suction  side  of  the  adjacent 
blade.  At  open  throttle  setting  a  strong  normal  shock  was  present  at  the  trailing  edge  of 
the  blade.  This  was  contrasted  by  the  near-stall  throttle,  Figure  52b,  where  the  shock  had 
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detached  from  the  blade,  which  created  an  unpinned  bow  shock.  The  oblique  and  normal 
shocks  that  had  moved  forward  of  the  passage  now  experienced  a  strong  tip  leakage 
vortex  were  it  interacted  with  the  detached  bow  shock  and  impinged  on  the  pressure  side 
of  the  adjacent  blade. 

Figure  53  shows  the  progression  of  the  casing  pressure  field  in  the  rotor  at  100% 
speed  in  the  rotor-only  configuration  as  it  approached  stall,  Block  A,  from  an  initial  full 
open  throttle,  Block  G. 


Figure  53.  Rotor  Only  Steady  State  Pressure  Contour  Plots  at  100%  Rotor  Speed 


Using  Figure  52  as  a  key  the  following  observations  of  Figure  53  can  be  noted: 

•  With  the  throttle  fully  open,  Block  G,  presented  a  oblique  and  attached 
shocks  with  relatively  weak  pressure  gradients.  At  this  throttle  setting 
nearly  all  the  compression  was  accomplished  by  the  flow  turning  through 
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blades  and  the  strong  normal  shock  at  the  trailing  edge;  this  corresponded 
to  the  low  efficiency  of  the  rotor  in  this  operation. 

•  Block  E,  displays  a  weak  attached,  a  stronger  oblique  shocks,  and  a  weak 
normal  shock  at  the  trailing  edge.  With  the  compression  being  shared 
most  effectively  between  the  obliques,  normal  shock,  and  blade  this 
corresponded  to  the  peak  effieciency  achieved  in  this  test  series. 

•  Block  D,  shows  a  sudden  disappearance  of  the  oblique  and  normal  shocks. 
This  marked  the  point  were  the  compressor  efficiency  has  begun  to  fall. 
This  was  also  the  point  were  the  tip-leakage  vortex  first  appeared  high 
velocity  air  was  injected  deep  into  the  blade  passage  which  deformed  the 
shock. 

•  Block  A,  the  pre-stall  condition  of  this  rotor  displayed  no  oblique  and  only 
a  normal  bow  shocks.  Also  present  was  a  strenghtened  tip  leakage  vortex 
that  impinged  on  the  pressure  side  of  the  blade. 

Figure  54  shows  the  progression  of  the  rotor  in  the  full  stage  configuration  at  95% 
rotor  speed  as  it  approached  pre-stall  conditions,  Block  A,  from  an  initial  full  open 
throttle,  Block  I.  The  Block  lettering  is  irregular  in  this  figure  so  that  labelling  would 
mirror  the  additional  figures  found  in  the  appendices.  It  is  also  important  to  note  that 
these  plots  differ  from  those  shown  above  with  the  elimination  of  the  Kulite  Probes 
positioned  at  the  trailing  edge  of  the  blade.  This  has  the  net  affect  of  creating  a 
“bleeding”  of  the  pressure  contours  after  the  last  probe. 


57 


Figure  54.  Full  Stage  Steady  State  Pressure  Contour  Plots  at  95%  Rotor  Speed 


While  Figure  54  maintains  a  similar  progression  from  the  attached  shock  and 
oblique,  Blocks  D  thru  I,  to  the  bow  shocks  seen  in  Blocks  A  thru  C  there  are  differences 
in  the  last  two  pre-stall  datasets: 

•  In  Blocks  A  and  B  the  characteristic  tip  vortex  had  diminished  vice 
strengthened  like  the  rotor-only  configuration.  In  its  place  was  a  more 
evident  shock-boundary  layer  interaction. 

This  difference  was  most  likely  due  to  the  flow  restriction  caused  by  the  stator 
that  had  the  effect  of  moving  higher  pressure  air  forward  into  the  rotor  blade  passages.  It 
is  important  to  note  that  a  similar  progression  from  open  throttle  to  pre-stall  occurred  at 
all  rotor  speeds  in  both  the  rotor-only  and  stage  configuration. 
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D.  CORRELATION  OF  FREQUENCY  AND  TIME  DOMAIN  ANALYSIS 

Drawing  together  the  products  of  the  time  and  frequency  analysis,  in  the  steady 
state,  an  attempt  to  correlate  the  findings  was  critical  in  providing  some  justification  for 
the  large  abrupt  changes  in  flow  that  were  recorded. 


Figure  55  shows  blade  passage  plots  for  the  TCR  at  100%  Speed  in  the  Rotor- 
Only  configuration. 
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Using  components  from  the  previously  discussed  analysis,  the  following 
observations  can  now  be  made: 

•  An  increase  in  RMS  power  from  Condition  C  to  Condition  B  was 
experienced  by  the  rotor.  This  corresponded  to  the  speed-throttle  profile 
where  the  bow  shock  preceding  the  blades  suddenly  detached  from  the 
leading  edge. 

•  A  decreased  in  RMS  power  from  Condition  B  to  Condition  A  was  also 
shown.  As  shown  in  Figure  55,  there  was  a  corresponding  forward  shift 
of  the  bow  shock  and  a  reduction  in  the  pressure  gradients  associated 
with  the  tip  leakage  vortex. 

These  observations  are  reinforced  by  Figure  56  where  the  standard  deviation  of 
the  pressure  was  plotted  within  the  blade  passages  for  the  same  throttle  settings.  As 
shown  there  was  a  characteristic  shift  of  the  high  deviation  oscillations  from  the  blade  tip 
to  mid-passage.  This  indicated  movement  in  the  bow  shock,  which  was  not  unusual,  as  it 
was  detached  from  the  leading  edge. 
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Similarly,  Figure  57  shows  the  progression  of  conditions  within  the  blade  passage 
at  95%  Speed  in  the  Full  Stage  Condition. 


Flow  Rate 


Figure  57.  Full  Stage  Pressure  Contour  and  Frequency  Data  Comparison  at  95% 

Rotor  Speed 


The  following  observations  are  made: 

•  An  increase  in  low  frequency  oscillations  immediately  after  the  attached 
shock  detached  from  the  blade  (Condition  D  to  C). 

•  An  attenuation  of  the  oscillation  after  the  higher  pressure  fluid  pushed 

the  shock  forward  ahead  of  the  blade  (Condition  C  to  D). 
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Figure  58.  Full  Stage  Standard  Deviation  Pressure  Contours  at  95%  Rotor  Speed 

The  standard  deviation  plot  for  these  runs,  Figure  58  shows  the  same 
characteristic  forward  shift  of  the  shock  and  reduction  in  tip  vortex. 

E.  STALL  INVESTIGATION 

As  discussed  previously,  stall  events  are  investigated  in  a  single  test  run  at  a  given 
speed  and  usually  requires  a  longer  time  run  out  than  the  standard  2-second  sample 
required  for  steady  state  operations.  Because  of  the  extreme  file  sizes  associated  with 
these  stall  tests  a  targeted  approach  was  required  so  that  only  those  rotations  of  interest 
were  analyzed. 

1.  Rotor  Only  Configuration 

Since  stall  was  always  accompanied  by  large  changes  in  rotor  speed  this  data  was 
plotted  so  that  the  rotation  where  stall  first  occurred  could  be  easily  identified. 
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Figure  59  shows  the  rotor  in  full  stage  configuration  as  it  approached  stall  at  95% 
rotor  speed.  As  shown,  there  was  a  characteristic  surge  in  rotor  speed  at  the  instant  of 
stall.  Since  the  initiation  of  a  stall  cell  was  associated  with  great  deviation  in  flow  as  it 
was  forced  around  the  blockage  into  adjacent  bade  passages  the  bulk  inlet  angle  was  then 
plotted  to  assist  in  the  identification  of  the  precise  rotation  where  stall  first  occurred, 
Figure  60. 


Figure  60.  Full  Stage  Bulk  Inlet  Angle  Entering  Stall  Event  at  95%  Rotor  Speed 
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Pressure  contour  plots  of  the  entire  rotor  were  then  assembled  over  the  identified 
rotations  in  the  Stead,  Pre-Stall,  and  Stall  regions  identified  so  that  the  onset  conditions  of 
stall  could  be  more  clearly  visualized,  Figure  61.  The  features  that  characterize  these 
regions  within  a  full  stage  are: 

•  Blade  Passage  Instabilities  existed  in  all  operating  conditions;  including 
steady  operating  conditions. 

•  As  the  rotor  enters  stall  it  went  through  several  cycles  of  stall  cell 
formation  and  recovery.  While  the  stall  cell  grew  across  more  blade 
passages  it  moved  around  the  annulus  backward  in  the  relative  frame. 

•  Once  the  entire  rotor  was  stalled,  it  began  to  surge  where  it  partially 
recovered  and  stalled  as  discussed  in  the  introductory  chapter. 


Figure  61 .  Full  Stage  Entering  Stall  at  95%  Rotor  Speed 
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An  analysis  of  the  blade  passages  in  steady  state  operations  focused  on  the  region 
of  instability  and  adjacent  effected  zones,  Figure  62.  As  shown  regions  adjacent  to  the 
instabilities,  that  decreased  incident  angles,  tended  to  “grow”  larger  tip  leakage  vortices. 
This  was  expected  since  this  change  in  incident  worked  to  further  decrease  pressure  on 
the  suction  side  of  the  blade  and  increases  the  gradient  across  the  blade  tip  driving  more 
flow  toward  the  blade  edge. 


Figure  62.  Full  Stage  Steady  State  Blade  Passage  Instabilities  at  95%  Rotor  Speed 


2.  Full  Stage  Configuration 

The  above  was  analysis  conducted  at  100%  of  rotor  speed  in  the  rotor  only 
configuration.  Plotting  the  rotor  speed  versus  time  revealed  the  approximate  location  of 
the  initiation  of  stall  at  the  instant  of  rotor  over-speed,  Figure  63. 
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This  approximate  location  was  then  used  to  identify  the  exact  moment  of  stall 
using  the  bulk  flow  inlet  angle. 


Stall 


Figure  64.  Rotor  Only  Bulk  Inlet  Angle  Entering  Stall  Event  at  100%  Rotor  Speed 


As  shown  in  Figure  64,  the  same  three  regions  of  steady,  pre-stall,  and  stall  were 
still  present  however  the  region  of  pre-stall  was  greatly  diminished  in  length  and  did  not 
display  the  same  prominent  oscillations  that  the  full  stage  did. 
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Generating  contour  plots  of  rotations  within  three  regions  yielded  the  following 
characteristic  blade  passages  over  time,  Figure  65. 


Figure  65.  Rotor  Only  Entering  Stall  at  100%  Rotor  Speed 


An  analysis  of  the  rotations  within  these  operating  regions  yielded  the  following 
observations  when  compared  to  the  full  stage  analysis: 

•  The  rotor  still  exhibited  blade  passage  instabilities  even  when  operating 
in  steady  state  conditions. 

•  Inception  of  stall  was  initially  focused  around  a  few  blades  but  quickly 
began  to  expand  to  adjacent  blades  as  flow  was  redirected  around  the 
stall  cell.  The  relative  motion  of  the  stall  cell  was  retrograde  to  rotor 
motion. 

•  Once  the  stall  cell  expanded  to  cover  the  entire  blade,  it  did  not  go 
through  a  cycle  of  recovery  and  stall  like  the  full  stage. 
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An  analysis  of  stall  cell  formation  and  growth  in  the  steady  operations  region  was 
also  conducted,  Figure  66. 


Figure  66.  Rotor  Only  Steady  State  Blade  Passage  Instabilities  at  100%  Rotor  Speed 

Figure  66  shows  the  location  of  steady  state  instability  and  the  adjacent  blade 
passages  that  are  affected  by  the  flow  moving  around  the  blockage.  As  noted  in  the  full 
stage  instability  analysis,  areas  of  decreased  incident  tended  to  display  larger  tip  leakage 
vortices  because  of  the  low-pressure  regions  that  were  created  on  the  suction  side  of  the 
blade.  This  was  contrasted  by  the  areas  of  increased  incident  were  lower-pressure 
regions  were  formed  on  the  pressure  side  of  the  blade,  consequently  reducing  the  size  of 
the  tip  leakage  vortices. 
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VIII.  CONCLUSIONS 


The  purpose  of  this  research  was  to  develop  an  integrated  easy  to  use  data 
analysis  tool  to  enable  rapid  and  consistent  data  analysis  regardless  of  configuration  or 
rotor  geometry.  This  tool  then  enabled  a  detailed  investigation  of  the  Transonic 
Compressor  Rig  (TCR)  to  characterize  differences  in  the  behavior  of  the  rotor-only  and 
full-stage  configurations,  to  identify  candidate  pre-stall  indicators  that  could  be  used  to 
predict  stall  prior  to  its  inception,  and  to  develop  a  more  complete  understanding  of  how 
a  rotor  departs  from  the  idealize  axi-symmetric  assumption  when  operating  at  steady- 
state. 

A.  INTEGRATED  DATA  ANALYSIS  TOOL 

The  TCR  Data  Analysis  tool  developed  for  this  research  enabled  a  comprehensive 
analysis  of  all  the  data  captured  during  the  experimental  phase.  The  capability  to  rapidly 
access  the  relevant  information  within  a  given  test  within  hours  of  its  completion  will 
allow  future  experimentation  to  be  more  responsive  to  failures  in  equipment  or  enable 
deeper  investigation  of  unexpected  behavior  on  the  fly.  This  code  base  was  developed  in 
such  a  way  that  data  from  any  rotor  can  be  quickly  imported  and  analyzed  with  small 
adjustments  to  geometric  constants  available  in  the  GUI.  This  provides  a  consistency  to 
future  analysis  that  up  to  this  point  was  not  possible. 

B.  STEADY-STATE  OPERATION 

The  analysis  of  the  stage  in  steady-state  operation  revealed  a  configuration  that 
was  much  stiffer  than  the  rotor  only  configuration  since  the  RMS  power  was  less 
concentrated  in  the  low  frequency  bands.  This  stiffening  is  largely  explained  by  the 
stator  acting  to  effectively  reduce  the  air  mass  behind  the  rotor  since  it  was  this  mass  that 
is  inversely  proportional  to  frequency  which  is  raised  and  then  attenuated. 

A  strong  correlation  between  bulk  flow  maximum  relative  inlet  angle  and  stall 
inception  was  also  found  with  both  configurations  achieving  a  71 -degree  inlet  flow  angle 
in  the  sonic  speeds  (90%  to  100%)  and  75-degrees  inlet  flow  angle  in  the  subsonic  and 
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transonic  speeds  (70%  and  80%).  Further  analysis  comparing  the  frequency  domain  data 
with  the  time-average  pressure  contour  plots  identified  a  connection  between  RMS  power 
and  the  presence  of  an  unpinned  bow  shock. 

Investigation  of  the  standard  deviation  measured  at  each  data  point  within  the 
passage  over  an  entire  test  run  also  showed  that  regions  of  large  deviation  traveled 
forward  of  the  blade  and  into  a  mid  pitch  location  just  prior  to  stall.  In  all  cases,  this 
instability  travelled  from  the  blade  tip  toward  0.6  times  the  blade  pitch.  These  two 
factors,  bulk  inlet  angle  and  instabilities  moving  toward  pitch  center,  stood  out  as 
indications  that  a  stall  event  was  imminent. 

C.  STALL  INCEPTION  ANALYSIS 

Analysis  of  the  rotations  just  preceding  and  following  a  stall  event  focused  on  the 
steady  nature  of  any  pre-stall  instability,  stall  cell  development  and  growth,  and  the 
behavior  of  the  cells  once  the  cell  has  formed  across  the  entire  rotor.  In  both 
configurations,  the  stall  cells  formed  in  localized  region,  across  3-4  blade  passages, 
which  rotated  at  part  speed.  Once  the  throttle  position  was  adjusted,  the  stall  cell  grew  to 
cover  the  whole  annulus  within  5-10  rotations. 

Primary  differences  in  behavior  between  the  rotor-only  and  full  stage 
configuration  was  in  the  transition  period  from  steady  operation  toward  stall  and  the  stall 
event  itself.  The  stage  exhibited  strong  oscillations  in  the  inlet  angle  going  into  stall  and 
displayed  a  greater  ability  to  recover.  Once  having  stalled  it  then  began  to  partially 
recover  and  stall  in  a  cyclical  way.  The  rotor  by  contrast  showed  no  oscillation  in  inlet 
angle  during  the  transition  phase  and  displayed  no  ability  to  recover  from  stall  once  the 
annulus  was  completely  stalled. 
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IX.  RECOMMENDATIONS 


As  discussed  in  previous  theses  the  incorporation  of  a  more  realistic  engine 
installation  would  lend  more  validity  to  the  results  discussed  in  this  research.  This  would 
require  a  shorter  or  serpentine  intake  to  the  TCR  so  that  it  would  operate  similarly  to  any 
future  military  aircraft. 

Future  tests  using  the  full  stage  will  also  require  an  increase  in  the  number  of 
pressure  transducers  over  the  stator  so  the  interaction  between  the  rotor  and  stator  can  be 
better  evaluated  and  understood.  As  currently  outfitted,  it  is  not  possible  to  sense  any 
flow  disturbance  downstream  of  the  rotor  or  the  presence  of  instabilities  on  the  stator. 
Since  it  is  these  rotating  instabilities  on  the  stator  that  are  the  root  cause  of  compressor 
surge  identifying  their  location  and  strength  in  relation  to  stall  cells  on  the  rotor  is  critical 
to  developing  a  better  understanding  of  the  TCR  in  this  configuration.  This  system 
should  also  sample  up  to  the  Kulites  theoretical  limits  of  300  kHz  providing  a  greater 
amount  of  resolution  to  the  data  already  being  processed. 
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APPENDIX  A:  ROTOR-ONLY  KULITE  DATA 
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Standard  Deviation  Plots 
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APPENDIX  B:  FULL  STAGE  KULITE  DATA 


Full  Stage  Data 
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APPENDIX  C:  MATLAB  INTERFACE 


To  better  integrate  the  analysis  routines  used  in  this  research  a  MATLAB  GUI 
was  developed  using  the  Guide  GUI  Builder  so  that  new  datasets  could  be  quickly 
initialized  and  analyzed. 

Overview  of  GUI  Components 

As  shown  below  the  TCR  Data  Analyzer  GUI  is  separated  into  6  sub  windows. 
These  sections  control  specific  analysis  variables  (Window  1),  Environmental  and 
Equipment  Constants  (Window  2),  Test  Run  Constants  (Window  3),  Subroutine  Controls 
and  Selectors  (Windows  4  &  5),  and  analysis  initiation  and  interface  controls  (Window 
6).  Each  of  these  sections  will  be  described  below  in  the  context  of  a  full  analysis  cycle. 


Setting  up  MATLAB  Interface 

To  setup  your  specific  MATLAB  working  environment  to  use  the  TCR  Data 
Analyzer  GUI  you  must  first  copy  the  m-files  and  GUI  fig  file  to  the  directory  of  your 
choice  and  then  add  these  directories  to  MATLABs  path  list. 
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To  do  this  you  should  complete  the  following  tasks: 

1 .  Navigate  to  the  Set  Path  Menu  item  under  the  File  Menu. 

2.  In  the  new  window  you  must  choose  to  Add  with  Subfolders... 


3.  Navigate  to  the  directory  that  you  copied  the  GUI  files  to  and  select  OK. 

4.  Save  the  changes  and  Close  the  Set  Path  window. 

5.  Start  the  program  by  typing  TCR_Data_Analyzer  in  the  MATLAB  Command 
Window. 

Initializing  a  Dataset 

Prior  to  initiating  a  data  analysis  routine  the  GUI  must  first  be  given  the  test  run 
specific  parameters  such  as  Run  Numbers,  Kulite  Placement,  and  Test  Run  Title.  This  is 
accomplished  through  the  use  of  the  Test  Run  Constants  window,  below. 

To  use  this  interface  you  must  first  use  the  Select  Data  Set  Directory  Button 
from  Window  6  and  navigate  to  the  directory  where  the  Test  Run  CSV  files  are 
located. 

A  description  of  the  input  parameters  is  included  below: 
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Test  Run  Constants 
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•  Problem  Title  -  This  is  a  matter  of  convenience.  Usually  a  title  includes  the  rotor 
speed  and  configuration  although  notes  regarding  peak  efficiency  or  throttle 
setting  may  also  be  included. 


An  example  Title  for  a  100%  speed  run  in  the  Rotor  only  configuration  might  be 
“95%  Rotor  Only.” 
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Run  Numbers  -  These  values  come  from  the  sequence  of  the  test  run  and  usually 
consist  of  four  numbers  although  some  early  runs  had  three.  The  exact  run 
numbers  can  be  located  in  a  text  fde  in  the  root  of  the  test  run  that  is  named  after 
the  date  of  the  testing. 

Impact  Kulite  Probe  -  This  value  is  always  “6.”  The  ability  to  choose  a  different 
Channel  Number  is  to  allow  easier  upgrade  of  equipment  in  future  testing. 

Kulite  Probe  Data  -  This  table  contains  details  regarding  the  number  of  Kulites 
(For  plotting  purposes)  the  axial  and  radial  location  of  the  Kulite  Probes  and 
identifies  the  Static  Pressure  probe  that  should  be  used  for  calibration.  Default 
values  for  this  window  are  automatically  populated  although  changes  can  be  made 
as  needed  to  accommodate  changes  in  sensors  for  a  particular  series  of  test  runs. 
This  window  contains  the  functionality  to  tell  MATLAB  what  channels  are  being 
recorded  and  if  any  of  the  recorded  channels  should  be  excluded  from  the 
analysis.  As  shown  below,  a  sample  test  run  is  recording  channels  1  thru  6  and  10 
to  1 1.  This  window  also  shows  that  Channel  10  is  a  dead  channel  and  should  be 
excluded  from  the  analysis. 
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If  mistakes  are  made  in  entering  the  data  the  tables  can  be  reset  using  the  Reset 
Tables  Button  located  at  the  bottom  comer  of  the  window. 


•  Once  all  changes  have  been  made  you  should  then  Save  the  changes  to  the  dataset 
directory.  While  saving  changes  in  not  strictly  required  to  run  the  GUI  it  allows 
MATLAB  to  call  up  this  file  and  read  in  the  values  directly  if  further  analysis  is 
required  on  a  particular  dataset. 

Data  Analysis 

A  steady  state  analysis  consists  of  Frequency  and  time  domain  analysis  of  the 
Kulite  and  Hotwire  signals.  To  initiate  a  steady  state  analysis  the  following  procedure 
should  be  used: 

1.  The  dataset  should  be  initialized  using  the  initialization  procedure  described 
above. 

2.  Select  the  Analysis  from  the  Analysis  Type  Window. 


-  Analysis  Type- 

0  Steady  State  (a »  Steady  Contour 

Correct  Shock 

©  Freq  Contour  Transient  Contour 

□  HQ  Movie 


There  are  4  primary  analysis  methods  available  and  two  analysis  modifiers  that 
only  operate  on  their  parent  analysis.  These  4  modes  are: 

•  Steady  State  -  Includes  a  series  of  time  and  frequency  analysis  that  include 
output  plots  of  pressure,  inlet  bulk  properties,  and  water  fall  plots  of  Kulite 
Signal  FFTs. 

•  Steady  Contour  -  Produces  a  time  averaged  plot  of  a  single  passage.  Can  be 
modified  with  the  Correct  Shock  switch  in  the  event  that  the  shock  angle  is 
not  calculated  properly. 

•  Transient  Contour  -  Produces  a  pressure  contour  plot  over  an  entire  rotation, 
or  series  of  rotations,  similar  to  above  but  without  the  time  averaging.  Also 
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produces  a  movie  of  the  pressure  contour  over  time.  Can  be  modified  with  the 
HQ  Movie  switch  to  ensure  video  is  saved  without  compression  (Useful  for 
presentation  purposes). 

•  Frequency  Contour  -  Provides  a  time-averaged  plot  similar  to  the  Steady 
Contour  options  but  also  includes  the  standard  deviation  and  frequency 
response  of  each  data  point. 

3.  Select  the  data  types  to  be  reevaluated. 


- Reload  Data - 

m  All  Data  □  Kulrte  Signal 

Steady  Contours  Q  Kulite  Calibration 


•  This  is  an  optional  choice  in  the  event  that  data  needs  to  be  reloaded  or 
recalibrated.  Usually  only  necessary  if  changes  in  the  code  is  changed  or 
bad  Kulites  are  discovered  and  need  to  be  eliminated  from  the  analysis. 

4.  Modify  any  of  the  variables  in  the  Analysis  Variables  window  as  needed.  The 
shown  values  are  the  default  and  usually  don’t  require  adjustment. 
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Analysis  Variables 


Rev  Start 

Rev  End 

Rev  No 

-1 

-1 

Rev  Samples 

Over  Samples 

FFT  Length 

1 

1 

IS 

Analysis  Region 
0  :  2 

Contour  Size 
150 


•  Rev  Start  and  Rev  End  -  Allows  a  targeted  analysis  of  a  specific  segment 
of  the  sample.  For  analysis  of  data  sets  that  don’t  include  stall  run  out  test 
this  should  be  left  alone.  The  range  of  revolution  should  be  approximately 
1000  for  best  results. 

•  FFT  Length  -  related  to  the  resolution  of  the  frequency  analysis  plots. 

•  Analysis  Region  -  Specifies  beginning  and  end  integer  locations  for  signal 
filtering  and  waterfall  plot  range. 

•  Contour  Size  -  Specifies  interpolation  size  for  all  pressure  contour  plots. 


5.  Adjust  any  equipment  or  environmental  variables  as  needed.  These  are  self- 
explanatory  and  should  not  require  any  adjustment  unless  there  are  geometry 
changes  to  equipment. 
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Environmental  £  EqjUfNnenrt  Constants 


Blade  Angle 
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Chord  J_ength 

20 

22 

D.8B824 

Gamma  Gas 

Grew  Const 

Rho  Hg 

1.4 

9.81 

13558 

Hub  Dioaneter 

11 

Casing  KulMe  Probes 

1 

7 

8 

6.  Initiate  Analysis  by  pressing  the  Run  Analysis  Button. 


7.  If  desired  all  the  generated  figures  can  be  saved  to  a  directory  of  your  choice  by 
selecting  the  Save  Figures  button. 
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